A Strongly Dipolar Bose-Einstein Condensate of Dysprosium 
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We report the Bose-Einstein condensation (BEC) of the most magnetic element, dysprosium. 
The Dy BEC is the first for an open /-shell lanthanide (rare-earth) element and is produced via 
forced evaporation in a crossed optical dipole trap loaded by an unusual, blue-detuned and spin- 
polarized narrow- line magneto-optical trap. Nearly pure condensates of 1.5x10^ ^^^Dy atoms form 
below T = 30 nK. We observe that stable BEC formation depends on the relative angle of a small 
polarizing magnetic field to the axis of the oblate trap, a property of trapped condensates only 
expected in the strongly dipolar regime. This regime was heretofore only attainable in Cr BECs via 
a Feshbach resonance accessed at high magnetic field. 

PACS numbers: 03.75.Hh, 67.85.Hj, 37.10.De, 51.60.+a, 75.80.+q 



The interplay between emergent electronic or spa- 
tial crystallinity, magnetism, and superfluidity is cen- 
tral to some of the most interesting materials of late, 
e.g., cuprate and iron-based superconductors, strontium 
ruthenates, and solid helium [1 . Quantum degener- 
ate gases possessing strong dipole-dipole interactions 
(DDI) [2] are thought to provide access to strongly cor- 
related quantum phases involving quantum magnetism, 
spontaneous spatial symmetry breaking, and exotic su- 
perfluidity. Supersolid [3] and quantum liquid crystal [4] 
phases may be accessible using strongly dipolar Bose and 
Fermi gases; with the realization of such quantum gases, 
we may shed light on characteristics of these phases dif- 
ficult to observe in their condensed matter settings [1]. 

We report the first realization of a strongly dipolar 
quantum gas at low field, a Bose-Einstein condensate 
(BEC) of ^^^Dy, which is an atom with unsurpassed 
dipole moment d = lO/iB, where /i^ is the Bohr magne- 
ton [5]. For comparison, Rb and Cr's moment equals l/i^ 
and respectively. The Dy BEC is the first quan- 

tum gas of a highly complex, open /-shell lanthanide, 
and opens a new frontier for exploring scattering behav- 
ior involving DDTinduced universality [6 , electrostatic 
anisotropy [71 [8], and an open /-shell submerged un- 
der closed s-shells. These complex atomic properties ev- 
idently do not prevent Bose-condensation, and explor- 
ing such complex collisional physics will greatly aid at- 
tempts to understand collisional behaviors of ultracold 
polar molecules, whose Bose-condensation has yet to be 
achieved. 

Unlike the ^^Cr dipolar BEC [9], we show that 
the ^^^Dy BEC reaches the strongly dipolar regime 
edd=MoM^^/127rfi^as>l without careful minimization of 
the 5-wave scattering length using a high-magnetic 
field Feshbach resonance [10] . The use of Feshbach reso- 
nances can lead to detrimental three-body loss [TTj and 
precludes the study of ultracold dipolar physics near zero 



field. Moreover, ultracold dysprosium suffers no chemical 
reactions like certain polar molecules [12]. 

Short- wavelength optical lattices confining Dy are ca- 
pable of exploring uncharted strongly correlated phases 
beyond the familiar Mott insulator at half-filling. 
Namely, density waves of various filling factors and lat- 
tice supersolids may now be accessible [2| ^ [13] without 
multilayer lattice enhancement [14 . Moreover, creation 
of unconventional and anisotropic superfluids |15j as well 
as explorations of ID strongly correlated gases [16] , trap 
instabilities [2 , spin textures and topological defects [17], 
roton modes [18 , and emergent structure in layered dipo- 
lar gases [19 are among the phenomena now within the 
realm of experimental possibility. 

The cooling procedure follows aspects of the Sr, Ca, 
and Yb BEC experiments [20], and we only provide es- 
sential experimental details before briefly discussing some 
of the ^^^Dy BEC properties. The BEC of ^^^^Dy, natural 
abundance 28.2%, is formed by forced evaporative cool- 
ing in a crossed optical dipole trap (cODT). The cODT 
is loaded with ultracold atoms from an unusual, blue- 
detuned narrow-line magneto-optical trap (MOT) [21] 
formed on the 1.8-kHz 741-nm line [22 , the demonstra- 
tion of which for Dy we report for the first time. This 
red- wavelength MOT (rMOT) is loaded by a repumper- 
less blue- wavelength MOT (bMOT) formed on the broad, 
421-nm transition. A Zeeman slower and transverse cool- 
ing stage, both at 421 nm, slow and collimate an atomic 
beam from a 1250 °C oven. Details regarding the bMOT 
are provided in Refs. [23 , and a subsequent publication 
will thoroughly describe the Dy rMOT and cODT. 

Figure [l] depicts the relevant Dy energy level dia- 
gram with cooling and trapping lasers as well as a tim- 
ing diagram for the experiment. The bMOT captures 
5x10^ 164 atoms, 70% of which, in steady-state, are 
~1 mK metastable atoms concurrently confined within 
the bMOT's magnetic quadrupole trap. We operate the 
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FIG. 1. (color online), (a) The relevant ^^^Dy level scheme. The nuclear spin is / = and the total electronic angular 
momentum isJ = L + 5' = 8 with L = 6 and S = 2. Experimental setup (b) and typical timing diagram (c). The trapezoid 
and diagonal line in the RF row represent the RF power ramp and the 0.7 MHz sweep about 5.65 MHz, respectively. 



bMOT at 1.4x10^ total atoms for this work. Of these, 
7x 10^ are captured in the rMOT — whose beams are over- 
lapped with the bMOT's — during the loading phase last- 
ing 3.5 s. Concurrently, metastable Dy decays to the 
ground state with a time constant of 2.3 s~^ [23] . 

The rMOT differs from most narrow-line MOTs, e.g., 
Ref. [24], in that the cooling laser is blue-detuned from 
resonance. For heavy, highly magnetic atoms possessing 
a narrow optical transition, a spin-polarized MOT may 
form below the center of the quadrupole field such that 
gravitational, magnetic, and optical forces balance each 
other. This technique, first employed for an 8-kHz-wide 
line in Er (/i = Tfis) [2Tj, creates an ultracold, dense 
gas at a position determined by the laser detuning such 
that the magnetic field Zeeman-shifts the otherwise red- 
detuned maximally weak- field-seeking state (MWS) to 
the blue of the laser frequency. The cloud shifts position 
in the magnetic gradient as a function of laser detuning, 
but the laser linewidth need not be as narrow as the 
transition, as often the case for narrow-line MOTs. 

The Er experiment captured 10-20% of the 10^ atoms 
loaded into the Er bMOT [25^ by fixing the rMOT 
laser intensity and detuning while ramping down the 
quadrupole field before recompression. We improve 
the capture efficiency of this technique by sweeping 
the laser's blue-detuning from A741 = +2.6 MHz to 
+0.6 MHz while reducing the laser intensity (quadrupole 
gradient) from / = 26 W/m^ {VB^ = 20 G/cm) to 4.8 
W/m^ (3.9 G/cm) [26]. We capture 50% of the total 
bMOT-cooled population by following this procedure. 
We hold this configuration for 1.2 s to cool to 12 /iK 
and spin-purify the gas in mj = 8, as determined by 
subsequent Stern-Gerlach (SG) measurements. While ul- 



tracold, the gas does not reach the line's 84-nK Doppler 
limit; similarly, the Er rMOT's 2 /iK temperature was 
substantially larger than its 190-nK Doppler limit [2T] . 

Reference [21] suggests that the seemingly hotter tem- 
perature arises due to momentum kicks from imperfectly 
extinguishing the quadrupole field before time-of-flight 
(TOF) measurements. However, these hotter tempera- 
tures may also arise from dipolar relaxation heating. In- 
deed, the plain evaporation of Dy after optical dipole 
trap loading implies the gas is no colder than 10 /iK 
when in the rMOT. Inelastic dipolar relaxation rates 
scale strongly with dipole moment and are expected to 
be quite rapid for Dy and Er [27] . We speculate that the 
gas heats by releasing Zeeman energy as the MWS state 
relaxes to lower-energy Zeeman states, while the rMOT 
lasers re-cool and repump to the MWS state. An equi- 
librium in temperature and spin polarization is reached 
for a given V^^, A741, and /, and future numerical sim- 
ulations will investigate this process [28 . 

We employ two horizontal ODTs, ODTl and 0DT2, 
to form the cODT at a position beyond the edge of 
the rMOT. Intersecting the ODTs with the rMOT dur- 
ing loading repels the gas due to a ~10-kHz Stark shift 
of the 741-nm line, which repositions the atoms in the 
quadrupole field to compensate the change in effective 
rMOT detuning. The optimal loading strategy thus 
involves setting the 18 W ODTl's propagation height 
roughly r = 1 mm below the rMOT center and fine- 
tuning the transfer efficiency with the 741-nm laser fre- 
quency while preserving spin-polarization and minimiz- 
ing gas temperature. 

The transfer efficiency from the rMOT to ODTl, as 
well as the initial gas temperature in the ODTl, is very 
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sensitive to the value of r. During the rMOT cool- 
ing stage, 20-kHz drifts of the 741-nm laser frequency 
(linewidth 20 kHz) shift the cloud in z enough to de- 
crease transfer efficiency. A Rb saturation-absorption 
reference laser stabilizes a transfer cavity to which the 
741-nm laser is locked; together this provides sufficient 
frequency stability for BEC production. 

Following this procedure, we obtain up to 4x 10^ atoms 
in ODTl at a temperature of 10 /iK, which is 15 x 
colder than the trap depth of ODTl. We experimen- 
tally determine the polarizability of Dy at 1064 nm to be 
116 a.u. (atomic units) by measuring the ODTl trap fre- 
quencies. This value is 65% smaller than recently calcu- 
lated [29 and will be used to improve magic wavelength 
estimates for future Dy optical lattices. 

ODTl is derived from a 30-W 1064-nm fiber laser, 
while 0DT2 is derived from a 25-W diode-pumped 
Nd:YV04 1064-nm laser. Beams from both are intensity- 
controlled by acousto-optic modulators (AOMs), and 
each beam is shaped into a cylindrical waist (e~^ radius) 
focused at the cODT center. The ODTl has a waist of 
30 (60) /im and trap frequency of 1 (0.5) kHz in z (p), 
where p lies in the x-y plane. 0DT2 has a waist of 22 
(70) /im and trap frequency of 1.2 (0.4) kHz in z (p). 

After extinguishing the rMOT lasers and quadrupole 
field, we immediately spin-flip the atoms from the MWS- 
state to the maximally strong-field seeking (MSS) ground 
state to prevent heating from dipolar relaxation. No op- 
tical pumping is needed because, as mentioned above, the 
rMOT— unlike a traditional MOT— provides a MWS- 
state spin-polarized gas [21]. RF adiabatic rapid pas- 
sage (ARP) flips the spin within 2 ms from mj = 8 to 
vfij = — 8 by using a = 0.3 G magnetic field and 
a swept-frequency RF-driven coil along z. Stern- Gerlach 
measurements indicate that ARP produces a nearly pure, 
absolute ground state gas in mj = — 8, and we use such 
ARP sweeps and SG measurements to calibrate the ap- 
plied magnetic fields and to zero the ambient magnetic 
field to <20 mG. The 0.3 G field (Am = 1 Zeeman 
energy equal to 25 /iK) is then rotated to point along 
z and increased to 5^ = 2.3 G (190 /iK) to suppress 
thermally driven spin transfer to mj > —8 states. 

Within 4 s of loading Dy from ODTl into the 340 /iK 
cODT, formed by ramping 0DT2 to 12 W in 1 s, plain 
evaporation both reduces the atom number to 1x10^ at 
the cODT center and decreases the temperature to 2 /iK. 
The trap frequencies are [fx, fy, fz] = [380, 500, 1570] Hz 
before the forced evaporation sequence begins. Power 
in the ODTs are reduced using the functional form Pf = 
Pi/ {l^t/rY [30 using experimentally optimized param- 
eters /3i = 2 and r = 15 s for ODTl. After a 9 s delay, 
0DT2's power is reduced for 23 s using j32 = 1.3 and 
r = 15 s. To recompress the trap and obtain an oblate 
shape — fz/fp > 1, with fp = y/ fxfy — both beams are 
linearly ramped in 0.8 s to their final powers. For the 
data presented here, final ODT powers yield the trap 
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FIG. 2. (color online). TOF profiles of the spin-purified Dy 
gas for three evaporation time- const ants, with r = 15 s in (e) 
and (f). (a,c,e) Data at centers are fit to a parabolic profile 
(green), which underestimates condensate fraction, whereas 
the distributions' wings are fit to a Gaussian profile (red). 
(b,d,f) Absorption images of the emerging BEC. (b) The 
transition temperature is 99(5) nK, with condensate fraction 
14.4(3)%; (b) 44(2)% condensate fraction at 56(3) nK; (c) 
a BEC of condensate fraction 73(4)% and 1.5(2) x 10^ atoms 
forms at 28(2) nK with density 10^^ cm"^. 



frequencies [fx^fy^fz] = [205,195,760] Hz, which pro- 
vide a ratio f^/fp = 3.8 [31 . The BEC is held for 0.2 s, 
the first 0.1 s of which the field may be rotated into 
the p-plane. Time-of-flight measurements are made by 
rapidly extinguishing the cODT and, after 1 ms, rotating 
the field into the imaging axis [see Fig. [ijb)]. Absorp- 
tion images are obtained using resonant 421-nm light and 
an exposure of 60 /is, much less than the decay time to 
metastable states [23 . 

The emergence of the Dy BEC below T = 100 nK is 
shown in Fig. [2] A dense parabolic peak emerges from 
a Gaussian background signaling the formation of a con- 
densate with high purity below T = 30 nK. The DDI 
preserves a parabolic condensate profile [32j , but the cri- 
terion for the Thomas-Fermi (TF) limit is modified for 
a dipolar BEC. Approximate expressions for the dipolar 
TF limit in the spherical and highly elongated prolate 
and oblate BEC regimes have been obtained [33 . Al- 
though our moderately oblate BEC fails to fall into these 
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regimes, we may use the oblate BEC analysis in Ref. [33] 
to estimate that the Dy BEC satisfies the dipolar TP 
limit for "effective" a^'s [34] in the vicinity of the recently 
estimated van der Waals length a = 76ao [8 . 

We calculate that an initial phase space density for 
BEC is reached at 120(20) nK for a trap with our mean 
frequency 310 Hz and atom number 7x10^. Applying 
corrections [9] for finite size effects and mean field en- 
ergy (using as = lOOao), we estimate a Tc = 100(20) nK, 
which is consistent with our measurement of 99(5) nK. 
The nearly pure BEC in panels (e) and (f) contains 
1.5x10^ atoms at a density of 10^^ cm~^. Further evi- 
dence of ^^^Dy's Bose-condensation is provided in Fig. [sj 
which shows the characteristic aspect ratio inversion of 
an interacting BEC released from an anisotropic trap. 
These data are obtained with a small, Bz = 2.3 G field 
parallel to the oblate trap's axis along z. BECs also form 
with 2.3 G fields rotated to an angle either 6> = 45° or 
60° with respect to the z-axis. However, no BEC forms 
in a 90° field {B^ = 0, Bp = 2.3 G), and a thermal cloud 
is obtained regardless of evaporation parameters. 

Strongly dipolar BECs — defined as Cdd > 1 — have 
been shown 21 [35H37 to be unstable for harmonic traps 
with 7 = f\\/f± < 1 but metastable for traps with 
7 = > where /y is the trap frequency parallel to 

the polarizing field, and f± is the geometric mean of the 
trap frequencies perpendicular to the field. The attrac- 
tive part of the DDI acts to elongate the cloud along the 
field direction, which if pointed along the trap's major 
axis leads to an instability much akin to that of non- 
dipolar BECs with large negative [2 . In contrast, 
aligning the field along the minor axis acts to suppress 
this elongation, minimizing the attractive contribution 
from the DDI to the mean field energy and stabilizing 
the gas. Weakly dipolar BECs — wherein Cdd < 1 — are 
stable regardless of the value of 7, because the repulsive 
contact interaction dominates the DDI [2 . 

While we observe a BEC in our oblate trap for ^ = 0° 
(7 = 3.8), the failure to observe a BEC in the oblate 
trap with 6 = 90° (7 = 0.51) is strongly suggestive of a 
BEC in the strongly dipolar regime [38]. Moreover, the 
tilt of the expanding condensate toward the polarizing 
field is indicative of a strongly dipolar condensate and 
was not observed in the weakly dipolar Cr system [39^. 
Gross-Pitaevskii equation simulations of the BEC stabil- 
ity boundary and tilt angle versus 6>, along with simula- 
tions of the TOF aspect ratio evolution, should provide 
a first measurement of ^^^Dy's effective [40 . 

In future work, we will search for broad Feshbach 
resonances and load the Dy BEC into optical lattices 
to search for phases predicted by the extended Bose- 
Hubbard model [2 . Quantum degenerate gases of Dy's 
other isotopes, including dipolar Bose- Fermi mixtures, 
could be formed with the experimental technique pre- 
sented here. Indeed, we recently optically dipole-trapped 
an ultracold dipolar Bose-Fermi mixture of ^^^Dy and 
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FIG. 3. (color online). Anisotropic expansion profile versus 
time after trap release, (a) and Vp are the dimensions of the 
parabolic profile fit to the BEC for ^ = 0°. Inset: Schematic of 
the oblate trap and magnetic field orientation, (b) Images of 
the expanding condensate after trap release. The condensate 
rotates by 7(1)° [9.4 (6)°] with respect to the ^ = 0° expansion 
orientation for = 45° [0 = 60°]. No BEC forms for 6> = 90°. 



^^^Dy — the first of its kind — using a modified version of 
this experimental technique; ongoing work strives to cool 
the mixture to quantum degeneracy. 
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